Time parameterization and stationary distributions in a relativistic gas 
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In this paper we consider the effect of different time parameterizations on the stationary velocity 
distribution function for a relativistic gas. We clarify the distinction between two such distributions, 
namely the Jiittner and the modified Jiittner distributions. Using a recently proposed model of a 
relativistic gas, we show that the obtained results for the proper-time averaging does not lead to 
modified Jiittner distribution (as recently conjectured), but introduces only a Lorentz factor 7 to 
the well-known Jiittner function which results from observer-time averaging. We obtain results for 
rest frame as well as moving frame in order to support our claim. 
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Introduction-Following a maximum entropy principle, 
Ferencz Jiittner [l| presented the first relativistic gen- 
eralization of MaxwcU-Boltzmann (MB) distribution in 
1911. Although Jiittner distribution has been gener- 
ally accepted and used in high energy and astrophysics 
0,11, 0, m, Ip, , some authors suggested several alter- 
natives |3.]lOr[lll [1^ that can be summarized in terms 
of the following r^-parameterized probability distribution 
functions (PDFs), 



(1) 



in which E — [m + p ) ' = 7717(7;) is the relativis 
tic single-particle energy with c 



Lorentz factor 7(1;) = (1 — v^Y^'^. 



1, rest mass m and 
Z is the normaliza- 
tion constant and /? is the temperature parameter. For 
?7 = and rj = \, the above PDFs reduce to the Jiittner 
function and the so called modified Jiittner function, re- 
spectively. 

The lack of rigorous microscopic derivations and ex- 
perimental evidences, made it difficult to decide which 
of the proposed relativistic distributions is the correct 
generalization of MB distribution. To resolve the uncer- 
tainty, semi-relativistic jl3| and fully relativistic [3, [lH| 
molecular dynamics simulations as well as Monte Carlo 
studies [1^ have been performed by different groups in 
recent years that unequivocally favored Jiittner distribu- 
tion. However, some recent investigations on relativistic 
Brownian motions [13, [3 have revealed that station- 
ary distributions can differ depending on the underlying 
time-parameterizations, a problem which never arises in 
Newtonian physics due to the existence of a universal 
time for all inertial observers. On the other hand, the 
maximum relative entropy principle (MREP) [l9l.[20l[2l| 
depicts how symmetry considerations lead to different 
stationary distributions, each with its own merits [2^ . 
Arguments of this kind suggest that one may possibly 



seek multiple relativistic generalizations of MB distri- 
bution. Using molecular dynamics simulations [23| and 
theoretical analysis [l3| , some authors have recently pro- 
posed that it is possible to establish a connection be- 
tween the time-parameter and the kind of symmetry that 
lead to a special distribution. Consequently, relativistic 
distributions fall into two classes that are distinguished 
by their associated symmetry and time-parameter. We, 
however, believe that the entire issue calls for a more 
careful reconsideration. 

In this article, we will first review how the choice of 
time parameters or reference measures affect the resulting 
stationary distribution. We further clarify the resulting 
consequence of these different choices. Using our recently 
proposed two dimensional hard disk model [l5| , we per- 
form simultaneous measurements at fixed observer's time 
'i' and particle's proper time 'r' with respect to a rest as 
well as a moving observer. The obtained results are then 
compared with stationary distributions consistent with 
different symmetry considerations, t-parametrization is 
shown to lead to Jiittner distribution, as is well-known. 
However, we find that r-parametrization, leads to a PDF 
that is consistent with a Jiittner function divided by a 7 
term, which is decidedly different from the original mod- 
ified Jiittner distribution obtained from MREP. This dis- 
tinction comes more to light when one considers a moving 
frame. 

Stationary distributions and symmetry considerations- 
The idea of maximizing relative entropy with respect 
to a pre-specified measure was a major step forward to 
clarify the underlying mathematical differences of the 
two mostly cited generalizations of MB distribution, i.e., 
Jiittner and modified Jiittner functions (2^ . 

Relative entropy [2l[ characterizes probability distri- 
butions with respect to a specific reference density. 
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in which, p(p) > plays the role of a reference density. 
It is then possible 22] to develop a general, coordinate 
invariant form of maximum entropy principle, under the 
constraints 
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d'p/(p) (3) 
e^Etot/N = J dV(p)i?(p). (4) 

The resulting stationary distribution is then given by, 

/(p)=p(p)exp(-/3i?(p))/Z, (5) 

Choosing a constant reference density, p{p) = po, the 
above equation reduces to the well known Jiittner func- 
tion. This choice correspond to a Lebesgue integrating 
measure, dp, = d'^p, which is the only translational in- 
variant measure in the momentum space [23|. Hence, it 
implies that Jiittner distribution is associated with a ref- 
erence measure that has translational symmetry. Since 
this kind of symmetry is not relevant in relativistic me- 
chanics, it is interpreted as a momentum conservation 
law during collisions. 

On the other hand, a momentum dependent reference 
density, p(p) = l/i?(p), leads to the modified Jiittner 
distribution. The integration measure consistent with 
this choice assigns to any subset A in momentum space, 
the measure number 



xiA) 



dp 



d'^p 



(6) 



Considering the fact that d'^p/E{p) is an invariant quan- 
tity under Lorcntz transformation, the modified Jiittner 
distribution is associated with a reference measure that 
is invariant under the fundamental symmetry group of 
relativity. 

It should be noted that we have no rigorous theoretical 
analysis or experimental evidence to favor one reference 
density /measure to the other. Numerical studies are also 
not decisive due to the arbitrariness in the choice of time 
parameters [13, HI] or discretization rules [23, [H, [2^ . 
Therefore, at this stage we are obliged to accept the ar- 
bitrariness in the choice of reference density. However, 
once the reference density is chosen, other parameters 
like Z and /? are deterministically obtained using the pre- 
specified constraints (Eqs. 3 and 2]). We will show how 
misleading conclusions can arise if one fails to appreciate 
this important point. 

Stationary distributions and time parameters- Dis- 
carding the classic notion of universal, absolute time in- 
troduces complexity to almost all time-dependent sub- 
jects in physics. Among these, are the evolution of dy- 
namical systems towards equilibrium and more specifi- 
cally their equilibrium probability distribution function 

[H, il, the descrip- 
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(PDF). As reported recently 
tion of the motion of relativistic particles with respect 
to the observer's time, t, and the particle's proper time, 
T are completely different. To elucidate, consider the 



t-averaged PDF of a Brownian particle in terms of the 
observer's time 



t 



dt'S{v-V{t')). 



(7) 



The velocity of the particle, V(i), can also be parame- 
terized in terms of the particle's proper time denoted by 
V(r). The corresponding r-averaged PDF is then defined 
as 



friv) 



1 



dT'5{v~V{T')). 



(8) 



By simply changing the integration variable, it can be 
shown [23 | that the relation between stationary {t, t — > 
oo) distributions is 



/oo(w) oc 



(9) 



The appearance of the factor 'y{v), in the denominator 
of Eq.® suggests that the r-averaged distribution is a 
modified Jiittner, if the t-averaged distribution turns out 
to be Jiittner function . Although the similarity 

is very suggestive, in this paper we show that a Jiittner 
function which is simply divided by E{p)(— mj(v)) is not 
necessarily equivalent to a modified Jiittner function. 

In order to distinguish the correct interpretation of 
Eq.®, we will check the validity of this equation by 
comparing the numerically obtained equilibrium distri- 
butions with respect to different time parameters, in the 
rest frame. Next, the results are further discussed by 
considering the problem from the viewpoint of a moving 
observer. In the light of these numerical simulations, one 
can decide whether the r-averaged distribution is truly 
described by a modified Jiittner function. 

The model we have used here is the previously pro- 
posed two dimensional relativistic hard-disk gas [l^. In 
this model, the disk-like particles move in straight lines 
at constant speed and change their momenta instanta- 
neously when they touch at distance a. The exchange 
of energy and momenta is governed by the relativistic 
energy-momentum conservation laws. In our simulation 
we have used N particles of equal rest masses m that 
are constrained to move in a square box of linear size L 
with periodic boundary condition. In order to simulate 
a stationary system in the rest frame, the center-of-mass 
momentum is put to zero manually. This condition would 
automatically be satisfied (if not at each instant but at 
least on average) if fixed reflecting walls were used [13] . 
The density is chosen to correspond to a dilute gas. 

In the next two sections we obtain f-averaged as well 
as r-averaged PDFs directly from simulations of such 
a model and compare them with the various proposed 
PDFs in order to clarify the relevance of various PDFs. 

Rest frame- To obtain the i-averaged PDF we let the 
system equilibrate (typically after 10^ collisions) and 
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FIG. 1; Equilibrium velocity distributions in the rest frame: 
numerically obtained x component of t-averaged (+) and r- 
averaged (•) velocity distribution from simulation of = 
100 particles of mass m = 0.1. Here, e = 3.06m and the 
corresponding temperature parameters are /3j = 7.62, /^ai j ~ 
4.86. A significant deviation from modified Jiittner function 
is evident. Similar results are obtained for f{vy). 



simultaneously measure velocities of all particles at a 
given instant of time t. To obtain the r-averaged PDF 
the proper-time of each particle, (i = 1, N), is com- 
puted during the simulation and velocities are recorded 
at a fixed proper-time value, = t. That is, veloci- 
ties are measured when the particles have the same life- 
time. To collect more data in the former (latter) scenario, 
we may either repeat the procedure for different initial 
conditions or perform measurements at several equally- 
separated time (proper time) instances. 

Simulation results for N = 100 particles are presented 
in Fig.l. Particles are initially placed on a square lattice 
of constant L/^/N and velocities are chosen randomly 
and the mean energy per particle is e = 3.06m. Velocities 
are all measured with respect to a rest observer at 2 x 10^ 
instants using time intervals lO^T, where T is the sys- 
tem's mean free time. As shown in Fig.l, the t-averaged 
(-I-) and T-averaged (•) PDFs agrees well, respectively, 
with Jiittner function with f3 — 7.62 (dotted red line) 
and Jiittner function with the same /3, divided by '£" 
(dashed blue line). This result confirms the relation de- 
scribed in equation Q. The more important question, 
however, is the correct interpretation of the right hand 
side of this equation. 

Considering the fact that the modified Jiittner func- 
tion differs from Jiittner function by a factor ^1/E\ one 
may conclude that r-averaged PDFs are described by 
modified Jiittner function in the same manner that t- 
averaged PDFs are depicted to be Jiittner function '23'|. 
Before drawing such a conclusion, we note that by substi- 
tuting modified Jiittner function in Eq.Q, one obtains a 
relation between the temperature parameter, /3, and the 
conserved quantities of the system (e,m), which is dif- 
ferent from the one obtained if one substitutes Jiittner 
function instead. For the two dimensional momentum 
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MJ 



(10) 
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where the indexes J and MJ refers to Jiittner and mod- 
ified Jiittner, respectively. In our simulation, this leads 
to Pj = 7.62 and /3mj — 4.86. Therefore, simply di- 
viding a Jiittner function by 7 does not lead to a modi- 
fied Jiittner function. The functional form as well as the 
temperature parameter are important in distinguishing 
a modified Jiittner function from its counterpart. This 
point is clearly demonstrated in Fig.l as the obtained 
PDF from r-averaging (•) does not fit well with a modi- 
fied Jiittner function (solid green line) but fits well with 
a Jiittner function divided by 7 (dashed blue line). 

Note that the difference of temperature parameters of 
Jiittner and modified Jiittner distributions in the above 
argument is a result of using same expectation values in 
the energy constraints in maximum relative entropy prin- 
ciple. It may be argued that this assumption is incorrect 
due to the fact that the two distributions refer to two 
different hyper surfaces in space-time [2^. Undoubtedly, 
a suitable choice of energy expectation value in Eq. pTj) 
gives the correct temperature parameter consistent with 
simulation data (i.e., /3 — (3j). The important question, 
however, is how should one obtain the new energy expec- 
tation value, consistent with Lorentz invariant measure, 
based on fundamental laws of relativity or statistical me- 
chanics? As far as maximum entropy principles are con- 
sidered, the relevant constraints should be specified with 
respect to the accessible information. Does changing the 
reference measure introduce new information that must 
be adapted as a constraint? 

At this stage we make no effort to answer these open 
questions and accept (at least with regard to numerical 
data) that the temperature parameter of both distribu- 
tion equals /3j. By this choice we still have two scenar- 
ios, first the r-averaged PDF is described by a modified 
Jiittner with P = f3j and second, it is a rescaled Jiittner 
function. These two are indistinguishable in the rest 
frame. However, because of different transformational 
properties of Jiittner and modified Jiittner functions un- 
der Lorentz boost, they take on distinctly different func- 
tional forms. We next consider our model in a moving 
frame in order to bias this distinction. 

Moving frame- To this end, we examine the system for 
an observer who is moving with a uniform velocity u in 
the negative x-direction with respect to the rest frame. 
Using the relative entropy maximization principle, the 
stationary distribution will be determined by an addi- 
tional constraint on the system, namely, that of a definite 
total momentum P' 1281. We therefore obtain: 
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exp[-/3j7(«)(£;'-u.p')] (12) 
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velocities separately. The dotted (red) and solid (green 
and pink) lines are the x{y) component of velocity dis- 
tribution obtained, respectively, by integrating Eqs. p^ 
and over Vyivx). We have also used the fact that 
the temperature parameters in the moving frame are the 
same as that of the rest frame [Illll. As clearly seen 
from numerical results, our two dimensional model shows 
that t-averaged PDF (-f ) is fitted to the expected x{y) 
component of Jiittner function (dotted red line) . The be- 
havior of T-averaged PDF (•), however, shows an evident 
deviation from the modified Jiittner function (solid green 
and pink lines). This result shows that no form of mod- 
ified Jiittner function can fit the numeric data obtained 
by proper time averaging. But how is a r-averaged PDF 
described in the moving frame? 

According to the arguments presented in the previous 
section, we consider a Jiittner function in the moving 
frame divided by E' , i.e., 



/'(pO - 



1 exp[-/3j7(»)(g^-u.p')] 
7(u)Z E' 



(14) 



FIG. 2: Equilibrium velocity distributions in the moving 
frame with relative velocity u — 0.1. The system param- 
eters are the same as Fig.l. In particular, /3j = 7.62 and 
/3a/j = 4.86. Part (a) shows the x component and (b) shows 
the y component distributions. The solid lines associated with 
modified Jiittner function with either /3j and (Smj do not fit 
the data well, while Jiittner function divided by E' with /3j 
fits perfectly well. 



rf _ 1 exp[- f3M.Jl{u){E' -u.p')] 

Jmj\P J — / x r, TTTcv r\ ■ y^"^ 

^{u)Zm.j 7(m)(£;' - u.p') 

The primed quantities are measured in the moving frame 
and the additional 7(m) term in the denominator, is 
due to the contraction of the moving box that encloses 
the system [2i|. Note that the reference density in the 
moving frame is obtained by p'(p') — Jp{v)i in which 
J = \dp^ /dp"^\ = det(a(^) is the Jacobian of the coordi- 
nate transformation, p'^ = a'^p'^ . Here, equals unity 
since for all proper Lorentz transformations, like boost 
and rotation, the determinant of the transformation co- 
efficients is det(a(^) = 1 [30]. 

Figure 2 shows simulation results for a system de- 
scribed in Fig.l with u = 0.1. Here, i-averaged PDF 
is obtained by measuring velocities simultaneously with 
respect to the moving observer as described in [3, [iBl ■ 
Considering the Lorentz invariance of proper-time it is 
not difficult to find the r-averaged PDF in the moving 
frame. In this case, the appropriate instant of measure- 
ment is the same as that of the rest frame, however, the 
velocities must be recorded as seen by the moving ob- 
server. To provide a better understanding of the effect of 
motion on PDFs we have shown the x{y) component of 



where Z is the normalization constant. As shown in 
Fig. 2, the x{y) component of the above distribution 
(dashed blue line), fits to simulation data as expected. 
These results again confirm that r-averaged PDFs are 
described by Jiittner functions which are simply divided 
by E. The similarity of such functional form to the mod- 
ified Jiittner function, especially in the rest frame, has 
recently led to misleading conclusions [l3, HI] . 

Concluding remarks- In the light of the above discus- 
sions one might now ask what is the correspondence be- 
tween symmetries that lead to different choice of refer- 
ence density in relative entropy with the symmetries that 
lead to the choice of time parameterization? In [2^ it is 
argued that the choice of coordinate time, i, and con- 
stant reference density results in a Jiittner distribution 
while a relativistically invariant measure along with a rel- 
ativistically invariant time (i.e., proper time r) leads to a 
modified Jiittner distribution. We, however, believe that 
the connection between various reference densities and 
time parameterization has not been well-established yet. 
Certainly a MREP should directly and uniquely lead to 
the equilibrium properties (i.e., PDF) once the reference 
density and constraints have been specified. To the best 
of our knowledge, there is no proof that choice of coordi- 
nate time leads to a constant reference density, just as no 
proof exists that choice of proper time leads to p = 1/E. 
Here, we have shown that if one calls what results from 
choosing different p's in the MREP as Jiittner and mod- 
ified Jiittner functions, the "modification" that results 
(due to the choice of p) in this method is different from 
the modification that results from a mere reparameteri- 
zation of time (which is in accordance with Eq.(9)). 

The authors kindly acknowledge the support of Shiraz 
University Research Council. 
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